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Gap solitons in optical lattices embedded into nonlocal media 
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We analyze the existence, stability, and mobility of gap solitons (GSs) in a periodic photonic 
structure built into a nonlocal self-defocusing medium. Counter-intuitively, the GSs are supported 
even by a highly nonlocal nonlinearity, which makes the system quasi-linear. Unlike local models, 
the variational approximation (VA) predicts the GSs in a good agreement with numerical findings, 
due to the suppression of undulating tails of the solitons. 

PACS numbers: 42.65.Tg, 42.65.Sf, 42.70.Qs 



I. INTRODUCTION 

Solitons are self-guided wave packets propagating in nonlinear media, maintaining the self-trapped shape. In 
particular, optical solitons are supported by the balance between the material nonlinearity and diffraction in the 
spatial domain or dispersion in the temporal domain [l|. As concerns spatial solitons in planar waveguides, it is 
well known that the self- focusing Kerr nonlinearity supports bright ones, while a defocusing nonlinearity admits 
dark solitons. Even in the absence of integr ability, solitons readily feature quasi-particle collisions, with the outcome 
depending on the relative phase between them. These properties suggest to use solitons in various applications to 
all-optical data-processing schemes and telecommunication systems [2|. 

Efficient control of the transmission and localization of light may be provided by photonic crystals (PhCs), built 
as structures with periodic modulation of the refractive index (RI). They open the ways to tailor the dispersion, 
diffraction, and routing of electromagnetic waves [3|]. Nonlinear PhCs, composed of appropriate materials, have 
revealed a wealth of nonlinear optical phenomena, including the self-trapping of localized modes in the form of the 
gap solitons (GSs) [4-16]. These modes can be formed in self- focusing and defocusing media alike, due to the possibility 
of the change in the sign of the effective dispersion/diffraction in PhCs. Q- Experimentally, GSs were first created 
in the temporal domain, as solitons in a short piece of a fiber Bragg grating "Sl. Technologies based on the use of 
reconfigurable (photoinduced) lattices, that have been implemented in photorefractive crystals [9[ and nematic liquid 
crystals [10|, offer new ways to control GSs in the spatial domain, by varying the lattice depth and spacings. 

Combining the benefits of PhCs and solitons. GSs have considerable potential for the use in photonics. GSs of matter 
waves have also been theoretically studied [ll| and experimentally created [i3| in Bose-Einstein condensates formed 
by atoms with repulsive interactions, trapped in optical-lattice potentials. Bifurcations and stability of optical GSs 
were analyzed in PhCs with the local Kerr nonlinearity ^IS] . However, the limited mobility of GSs in the transverse 
directions, due to their pin ning to the underlying lattice potentials [IJ], is an obstacle to the use of GSs in switching 
and routing operations |15l . [l6| . 

Recently, it has been predicted that solitons supported by a nonlocal nonlinearity, self-focusing or defocusing, in 
the combination with the effective diffraction induced by either the total internal reflection (ordinary solitons) [13] or 
bandgap spectrum [18l|, may move much easier across the lattice. Nonlocal effects come to play an important role as 
the characteristic correlation radius of the medium's response function becomes comparable to the transverse width 
of the wave packet jl9| . Experimental observations of nonlocal responses have been demonstrated in various media, 



including photorefractive crystals [20| , nematic liquid crystals [2lJ , and thermo-optical materials l22|.|23i . The nonlocal 
nonlinearity induces new features in the wave dynamics, modifying the underlying modulational [2J], azimuthal [25| . 
and transverse [26] instabilities. Suppression of the collapse of multidimensional solitons |27[, a change of interactions 
between them ,28,] , the formation of soliton bound states [2£] , merger of colliding solitons into a standing wave [30[ , 
and families of dark-bright soliton pairs |3l| were also predicted recently. 

The nonlocality is known to improve the stability of solitons due to the diffusion mechanism of the underlying 
nonlinearity. In the limit of strongly nonlocal nonlinearity, the system become an effectively linear one [32| . In such 
an extreme limit, the existence of GSs (for the defocusing sign of the nonlinearity) is questionable. In this work, we 
identify families of bright on-site and off-site GSs in self-defocusing nonlinear media by means of numerical methods 
and analytical methods. With the infinite range of the nonlocality, we demonstrate the existence of spatial GSs with 
a finite beam's width. The analytical consideration is based on the variational approximation (VA) with a Gaussian 
ansatz, similar to how it was applied to the matter- wave GSs in Refs. [3j, |3J|. Unlike the case of the local defocusing 
nonlinearity, in the nonlocal modes the Gaussian ansatz works well not only deep inside of the bandgap, but also close 
to its edge. The stability and mobility of the GS families in the nonlocal medium are investigated too. 

II. THE MODEL AND NUMERICAL RESULTS FOR GAP SOLITONS 

We consider a wave packet propagating along axis z in a PhC structure embedded into with a medium with the 
self-defocusing cubic nonlocal nonlinearity. A model widely adopted for the description of such media is [2J, [35| 

^^ = --^v, + y(,)v, + .,vi,, (1) 

n-d^n=\^\^, (2) 

where ^I* is the amplitude of the electromagnetic wave, x the transverse coordinate, n{x, z) a perturbation of the 
local RI corresponding to the intensity-response function with an exponential kernel, and d is a parameter which 
determines the degree of the nonlocality of the response. All the physical quantities and spacial coordinates are made 



dimensionless by normalization procedure with respect to the input beam width, wavelength, and Kerr coefficient 
of the nonlinear material [ij . The limit oi d ^ oo corresponds to the well-known Zakharov's system, which is a 
fundamental model in plasma physics (for Langmuir waves) and other fields [3^ . The PhC structure is represented by 
the periodic transverse potential, V{x) = Vosin x {x is normalized so as to make the period equal to tt). Stationary 
solutions with propagation constant /z are sought for as ^l/(x, z) = exp {—ifiz) 4>{x), which gives rise to the stationary 
version of Eqs. dH) and ([2|): 
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If the nonlinearity is omitted, Eq. ([3]) decouples from n{x) and becomes a linear equation, which supports Bloch-wave 
solutions, (^{x) ~ f (x) exp{i k x) , where k is the quasi-wavenumber, and f{x) is a periodic function with period tt. As 
an example, we take Vq = 4 and display the corresponding dispersion relation, including the three lowest three bands, 
in Fig. [1] From this diagram, it is seen that finite bandgaps are introduced by the periodic potential; in particular, 
the first finite bandgap covers a broad interval, 1.305 < /i < 3.19. 




FIG. 1. A typical example of the spectrum, with quasi-wavenumber k, induced by the linearize version of Eq. ([3]), with Vq = 4. 
Shaded areas are covered by the bandgaps. 



The nonlinearity may give rise to x-periodic modes [37[ or localized GSs [38|, with fi falling into the bandgaps, 
which refers to gap solitons. Starting with the GS solution in the middle of the gap, we have found different families of 
bright solitons numerically, using by standard relaxation technique with boundary conditions (f)(±oo) — 0. In Fig. [5J 
we demonstrate generic examples of the GS modes found in the first finite bandgap. Local (with d = 0), and nonlocal 
(for d = 40) on-site-centered GSs are shown near the bottom of the gap in Figs. [UJa, d), in the middle of the gap in 
Figs. [UJb, e), and close to the top edge in Figs. [2ljc, f), with propagation constants ^ = 1.31, ^ — 2.5, and fi = 3.1, 
respectively. Simplest higher-order GS solutions of the nonlocal model (off-site-centered solitons) are presented in 
Fig.Ufg-i). The two distinct types of the solitons, on-site and off-site, are defined by the position of their centers with 
respect to the underlying periodic potential [13, llMl • 

Similar to solitons in nonlocal media with the self- focusing nonlinearity [2J, [3^], the amplitude of the GSs in 
the present model increases at a higher degree of the nonlocality, d. As a result, the related total power, P = 
/_ l^(2;)P dx, is a growing function of d, as shown in Fig. |3l In comparison to the local nonlinear medium, with 
d = [ see Figs. [H^a-c)], the width of the RI perturbation, w„, becomes broader with the increase of d, for the focusing 
[2J,|35| and defocusing signs of nonlinearity alike, due to the diffusion type of the nonlocal response. Relations between 
Wn and d are shown in the first column of Fig. 2] On the contrary to the solitons in self-focusing nonlocal media 
|17lll8|. the beam's width of the GSs in the present case, Wb, decreases with the increase of d, as shown in the second 
column of Fig. 01 At very large values of d, the beam's width in the GS approaches a constant value. 
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FIG. 2. Typical examples of gap-soliton modes in the first finite bandgap. (a-c): On-site solutions in the local model (d — 0). 
(b-f): On-site solutions for d = 40. (h-i); Off-site solution for d — 40. Solid lines show the field profiles, while the corresponding 
profiles of the refractive-index perturbation, n{x), are plotted by dashed lines. 
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FIG. 3. The total power of on-site gap solitons versus the nonlocality parameter, d, for three distinct values of the propagation 
constant taken, respectively, near the bottom edge of the first finite bandgap's edge (/x = 1.31), in the middle of the gap 
(/I = 2.5), and approaching the top band edge (/i = 3.1). 



III. VARIATIONAL APPROXIMATION 



It is well known that undulating tails in the shape of GSs, induced by the underlying periodic potential, make the 
Gaussian ansatz inappropriate as an approximation for GSs, especially close to bandgap edges [33|. However, Fig. [2] 
demonstrates that the nonlocal nonlinearity makes GSs in the present model more localized, suggesting to apply the 
VA. The Lagrangian density for Eq. ^ is 
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FIG. 4. Left column (a,c,e): The width of the refractive- index profile, Wn, in gap solitons, versus nonlocality parameter d. 
Right column (b,d,f): The width of the field component of the gap solitons, w^,, versus d. Numerical and variational results are 
shown by solid and dashed lines, respectively. The value of the propagation constant is fixed in each panel, as indicated. 



Following the above argument, we adopt the Gaussian ansatz for field (\) and RI perturbation n, 



(a;, z) — j4(z)exp 



2wi{z) 



ib(z)x 



n{x, z) = C(z)exp 



2<(z) 



(6) 



with A{z\ b{z) and Wfc(z) standing for the amplitude, chirp and width of the field component of the GS, while C{z) 
and Wn{z) are the amplitude and width of its RI counterpart. Substituting the ansatz into Lagrangian density ([5]) 
and performing the standard calculations 14 1| . we arrive at VA-generated relations between the parameters, 



(3d + 2w2 



P (3d + 2wl) 2wn {3d + 2wl 



2wl 



8u;4 (2u)2 - dy 2wn {d + 2wiy 
2wl {2wl ~ d) 
3d + 2wl ■ 



{d + 2wlf 



(m + Vo), 



(7) 
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Using Eqs. (O and ([5]), in Fig. [S] we draw the surface plot for the width of the RI profile, w„, as a function of total 
power P and nonlocality parameter d, at a fixed value of /i, and compare it to the numerical results. As expected, 
the VA produces good results for the GSs taken in the middle of bandgap - for instance, at /i = 2.5. 

In addition, using the power and propagation constant found numerically in Sec. |TT1 we show in Fig. |4] that 
both numerical and variation solutions represent the same trend for the widths of both components of the GSs. In 
particular, for /i = 2.5 and /z = 3.1, as shown in Fig.|4l[c-f), Wn increases as the square root of the nonlocality strength, 
d, in agreement with Ref. |32| . On the contrary, Wb drops to a finite value as d increases. However, for jjl — 1.31, 
which is very close to the edge of the first finite bandgap, the trend is completely different. This difference is explained 
by the known fact that GSs with the propagation constant taken very close to edges of bandgaps are similar to the 
linear Bloch waves in the linear lattice, as seen in Fig. WiA)- When the amplitude of the respective undulating tails 
in the GS shape is comparable to its main peak, the Gaussian ansatz definitely fails. Nevertheless, close to the top 
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FIG. 5. (Color online) The surface plot for the width of the refractive- index component of the gap soliton, Wn, at /i = 2.5, as 
predicted by variational equations ([7}, ([§}. The chain of dots represents full numerical solutions. 



edge of the first finite bandgap - for instance, at /i = 3.1, the ansatz still works well because the major peak in the 
GS profile remains much higher than the undulating tails in the entire nonlocality regime. Thus, the applicability 
condition for the Gaussian-based VA in the system with the self-defocusing nonlocal nonlinearity is clear: It is usable 
as long as the GS propagation constant is taken not too close to the bottom edge of the first finite bandgap. 

In the extremely nonlocal regime, the field component in the GSs is much narrower than the RI profile. In this case, 
the RI profile may be approximated as n{x) « R{x) = e~^^^^^'^/{2\/d), where R{x) is the RI response function. The 
corresponding width of the RI profile is w„ ~ 2-\/dln2. Then, using the quasi-linear limit similar to that developed 
in Ref. [S^], one can predict the threshold power necessary for the formation of the GS, 



Pthr = 2(^ - ^lQ)Vd, 



(9) 



where fiQ as the propagation constant at the edge of the first finite bandgap. The comparison to the numerical results 
in this extremely nonlocal regime is shown in Fig. [B] 
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FIG. 6. The threshold power necessary for the formation of the gap solitons versus the nonlocality parameter, d. The solid 
and dashed lines depict, respectively, the numerical results and asymptotic approximation (|9]). 



IV. STABILITY AND MOBILITY OF THE GAP SOLITONS 

Having constructed the family of the GS solutions, we analyze their stability in the usual way, considering perturbed 
GS solutions as 



n = no + An, 



-iS* z'\ ibz 



(10) 

(11) 



where e ^ 1 is the perturbation amplitude , ua{x) is the unperturbed solution, and Im{(5} is the growth rate of the 
perturbations. Although the strengthening nonlocality makes the GS shape sharper and stronger localized, we have 
found, somewhat counter-intuitively, that the on-site GS family is still stable, while its off-site counterpart is not, cf. 
[13| . Figure [Tl^a) demonstrates that the nonlocality significantly reduces the growth rate of the unstable perturbation 
mode for off-site solitons, as in the case in the self- focusing nonlinearity [18|. Due to its diffusion character, the 
nonlocality smoothes down the undulating tails in the RI profile, n{x). It is this smoothness that stabilizes GSs, for 
either sign of the nonlinearity. Therefore, as the strength of the nonlocality increases, the GS solutions become more 
stable, through the broadening of the effective potential. 
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FIG. 7. (a) The instability growth rate and (b) the PN energy barrier versus the nonlocality strength, d, for different values of 
the total power, P. 



The mobility of the GSs in the present model was studied through calculating the respective Peierls-Nabarro (PN) 
potential barrier, which is defined as the height of an effective periodic potential generated by the underlying lattice. 
The potential barrier determines the minimum energy needed to move the center of mass of a localized wave packet by 
one lattice site [40]. As usual, we can calculate the PN barrier as the difference of values of the model's Hamiltonian 
(H) between on-site and off-site modes |13|, i.e.. 
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As seen in Fig. [Iljb), in the first finite bandgap the PN barrier is reduced in comparison to the case of the local 
nonlinearity, d = 0, which is a natural manifestation of the nonlocality. 



CONCLUSION 



We have reported the analysis of the existence, stability, and mobility of one-dimensional GSs (gap solitons) in 
the periodic potential structure combined with the self-defocusing nonlocal nonlinearity. We have found that the 
GSs become tighter localized in space, with a higher formation-power threshold. The results have been obtained 
in the numerical form and reproduced, with a reasonable accuracy, by the VA (variational approach). Using the 
linear-stability analysis and calculating the PN (Peierls-Nabarro) potential barrier, we have demonstrated that the 
GSs become not only more stable, but also more mobile, with the increase of the nonlocality. The comparison with 
the limit of the extreme nonlocality was reported too. Taking into regard the possibilities offered by the currently 



8 

available technology for fabricating nonlocal nonlinear media with controllable properties, such as photorefractive 
crystals, nematic liquid crystals and thermo-optical materials, the results reported in this work may suggest new 
possibilities for the design of soliton-based photonic devices. For instance, in photorefractive materials like SBN or 
LiNbOs, photorefractive gratings [r] and self-defocusing nonlinearity [42'| can be easily achieved by exteranl electrical 
field properly applied acorss the crystal axis. Also in liquid- filled photonic crystal fibers [4J|, 2D band gap effect with 
thermal-type self-defocusing nonlinearity is demonstrated. It may also be interesting to extend the model and the 
analysis of GSs in it (including vortex solitons) to the two-dimensional geometry. 
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